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Abstract 

We explicitly derive the principal Lyapunov Exponent in terms of the ra- 
dial equation of ISCO (Innermost Stable Circular Orbit) for Spherically symme- 
try(Schwarzschild, Reissner Nordstr0m) black-hole space-times. Using it, we show 
that the ISCO occurs at rjsco = 4M for extremal Reissner Nordstr0m black-hole 
and rjsco = 6M for Schwarzschild black-hole. We elucidate the connection be- 
tween Lyapunov Exponent and Geodesic Deviation Equation . We also compute the 
Kolraogorov-Sinai(KS) entropy which measures the rate of exponential divergence 
between two trajectories(geodesics)via Lyapunov Exponent. We further prove that 
ISCO is characterized by the greatest possible orbital period i.e. Tjsco > T p hoton 
among all types of circular geodesics(both time-like and null, geodesic and non- 
geodesic) as measured by the asymptotic observers. Therefore, ISCO provide the 
slowest way to circle the black hole. 



1 Introduction 

Nonlinearity of Einsteins equation makes nonlinearity of Einsteins General Theory of 
Relativity. So there may be a certain link between nonlinear Einsteins General Theory of 
Relativity and nonlinear dynamics. Particularly Lyapunov Exponent and Kolmogorov- 
Sinai(KS) entropy, are one of the bridges between them. In this paper, we shall focus on 
analytical calculations involving Lyapunov Exponents and Kolmogorov-Sinai in terms of 
equations of circular geodesies around a black-hole spacetime. This equatorial circular 
geodesies around a black-hole space-time playing a crucial role in General Relativity for 
classification of the orbits. It also determines important feature of the space-times and 
gives important information on the back ground geometry. 

The Lyapunov Exponent has been used to probe the instability of circular null 
geodesies and in terms of the Quasinormal modes for spherically symmetric space-time 
in [TJ, but the focus there is on null circular geodesies. It has been shown in this reference 
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that in the eikonal approximation, the real and imaginary parts of the Quasi Normal 
Modes of any spherically symmetric, asymptotically flat spacetime are given by(multiples 
of) the frequency and instability time scale of the unstable circular null geodesies. 

Note however that the Principal Lyapunov Exponents (A) have been computed in [5]- 
[7] using a coordinate time t where t is measured by asymptotic observers. Thus, these 
exponents are explicitly coordinate dependent, and therefore have a degree of unphys- 
icality. Here we compute the Principal Lyapunov Exponent (A) and Kolmogorov- Sinai 
entropy hk s analytically by using the proper time r which is coordinate invariant. Using 
them we study the stability of equatorial circular geodesies for Reissner Nordstr0m and 
Schwarzschild black-hole space-time. Another interesting point we have studied here is 
that the Lyapunov Exponent and Kolmogorov-Sinai entropy hk s can be expressed in terms 
of the ISCO equation and we find the ISCO occurs at r IS co — 4M for extremal Reissner 
Nordstr0m(RN) black-hole and rjsco — 6M for Schwarzschild black-hole via Lyapunov 
stability analysis and Kolmogorov-Sinai entropy hk s . 

The paper is organized as follows. In section 2 we provide the fundamentals of the 
Lyapunov Exponent. In section 3, the Lyapunov Exponent has been expressed in terms 
of the standard effective potential. In section 4 we elucidate the connection between Lya- 
punov Exponent and Kolmogorov-Sinai Entropy which is also determine the exponential 
divergence between neighbouring trajectories. In section 5 we derive the standard effec- 
tive potential in the equatorial plane. In section 6 we fully describe the equatorial circular 
geodesies, both time-like and null, for Reissner Nordstr0m blackhole and also discuss the 
Lyapunov Exponent can be expressed in terms of the ISCO equation and studied for the 
stability of the circular geodesies. In section 7 we discuss similar things for extremal 
Reissner Nordstr0m spacetime. In section 8 we relates the Geodesic Deviation equation 
in terms of Lyapunov Exponent for any spherically symmetric spacetime. In section 9 
we compare the ratio of angular velocity between ISCO and null circular geodesies. In 
section 10, we present our conclusions. 

2 Fundamentals of Lyapunov Exponent: 

The Lyapunov Exponents in a classical phase space are a measure of the average rates of 
expansion and contraction of a trajectories surrounding it. They are asymptotic quantities 
defined locally in state space, and describe the exponential rate at which a perturbation to 
a trajectory of a system grows or decays with time at a certain location in the state space. 
A positive Lyapunov Exponent indicates a divergence between two nearby geodesies, the 
paths of such a system are extremely sensitive to changes of the initial conditions. A 
negative Lyapunov Exponent implies a convergence between two nearby geodesies and 
the vanishing Lyapunov Exponent indicates the existence of periodic motion i.e. the 
system has a limit cycle. Lyapunov Exponents can also distinguish among fixed points, 
periodic motions, quasi-periodic motions, and chaotic motions. 
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An n-dimensional continuous-time(autonomous) smooth dynamical system governs by 
the differential equation [T3] of the form 

(1 T 

- = F(x;M). (1) 

The vector x consists n state variables, the function F describes the non-linear evolution 
of the dynamical system and M is a vector control parameter. The solutions are fixed 
point solutions or equilibrium solutions. Where t is defined as time parameter. The fixed 
points or critical points, are defined by the vanishing of the vector field; that is 

F(x;M) = 0. (2) 

Let the solution of (j2J) for M = M be x , where x G lZ n and Ai G 1Z m . To calculate 
the stability of this equilibrium solutions, we simply apply on it a small perturbation y 
and obtain 



x(t) = x + y{t) . 



(3) 



Substituting (j3J) into ([T]) yields 



dy 
dt 



F(x + y;M) . 



(4) 



Note that the fixed point x = Xn of (fll) has been transformed into the fixed point y = 
of (jl]). Expanding (J3J) in a Taylor series about xq and keeping only linear terms in the 
perturbation leads to 



J = F(x ;M ) + 



(5) 



or 



dx 



■y = Ay . 



(6) 



where the matrix A is called Jacobian matrix. If the components of F are 

x 2 , ar 3 , x n ), F 2 (xi, x 2 , x 3 , x n ), F 3 (x 1 ,x 2 , x 3 , ar n ) then 



/ 8F\ dF\ 
3ti dx'> 



A 



8x\ 



dF 2 9F2 
8x2 



8F n 8F n 
dxi 8x2 



8F\ x 

8Xn 



8F 2 
8x n 



(7) 



8F n I 

8x n ' 
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The eigen values of the constant matrix A provide information about the local stability 
of the fixed point x . The solution of the equation (jBJ) becomes 



y (t) = e^ A y . (8) 

where 

z — ' i\ 

i=0 

and yo e lZ n at time t G 7?.. If the eigen values A« of the matrix A are distinct, then 
there exists a matrix P such that P~ l AP = D, where D is a diagonal matrix with eigen 
values Ai, A2, ... ,X n ] that is 



D 



( Ai 
A 2 



\ 



\ 




A n / 



(10) 



If the eigen values are complex, then the matrix P will also be complex. The columns 
of the matrix P are the right eigen vectors pi,p2, ■■■,Pn of the matrix A corresponding to 
the eigen values Ai,A2,...A n ; that is, 



P 



fPi\ 

P2 



\Pn J 



(11) 



Hence 



AP 



( A Pl \ 

Ap 2 



\ A Pn J 



( XlPl \ 
A2P2 



\ KPn J 



PD 



(12) 
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Consequently 

D = P~ l AP . (13) 
Now transforming y = Pz into ([6]), we get 

P d 4 = APz 
dt 

or 

% = Dz. (14) 



dt 



Hence 



where 



z = e {t - to)D z . (15) 



z = z{t ) = P-'vo . (16) 

In terms of y, this solution becomes 

y{t) = PeWp- 1 ^ . (17) 

The matrix e^ - * ^ is a diagonal matrix entries e^* - * '**. Therefore the eigen values of A 
are also known as Characteristic Exponents or Lyapunov Exponents associated with F 
at (x ,M )- 

If we consider an initial deviation y(0), its evolution is described by 

y{t) = $(t)y(0) . (18) 

where $(t) is the fundamental(transition) matrix solution associated with the trajec- 
tory say X(t) which governs the dynamical equation ([1]). For an appropriate chosen y(0) 
in ( |18p the rate of exponential expansion or contraction in the direction of y(0) on the 
trajectory passing through X (trajectory at t = 0) is given by 



* - &G)Ki5ti) ■ (l9) 

where ||| denotes a vector norm. The asymptotic quantity A* is called the Lyapunov 
Exponent. 

If there exists a set of n Lyapunov Exponents associated with an n-dimensional au- 
tonomous system and they can be ordered by size that is 

Ai > A 2 > A 3 >,...,> A„ . (20) 
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The set of n-numbers Aj is called the Lyapunov Spectrum. 

Following Lyapunov|16j. the fundamental matrix $(t) is called regular if 



lim In | det$(t) | . (21) 



t— >oo 



exist and is finite and if there exists a normal basis of the n-dimensional state space such 
that 

n 

VA|= lim In | det9(t) | . (22) 

i=l 

If $(i) is regular, then according to a theorem Oseldec[2] the asymptotic quantity defined 
in (fl9|) exists and is finite for any initial deviation y(0) belonging to the n-dimensional 
space. 

The asymptotic quantity A', given by (1191) is also known as a one dimensional exponent. 
For p-dimensions, a p-dimensional Lyapunov exponent A is defined as 

'->°°W Vll!/i(0)A» 2 (0)A...A» p (0)||y 

where A is an exterior or vector cross product. 

The above fact of Lyapunov Exponents is valid for any non-relativistic system. But 
in General Relativity, there is no concept of absolute time, therefore the time parameter 
forces us to consider equation (TjQ) under spacetime diffeomorphism: u = ip(x), dr = 
rj(x)dt. As a result, chaos is a property of the physical system and does not depend 
on the coordinates used to describe the system. Therefore motivated by the work of 
MotterpQ chaos is characterized by positive Lyapunov Exponents. Similarly chaos is also 
characterized by positive Kolmogorov -Sinai Entropy (hl s ). They are coordinate invariants 
and transform according to 

< V >t 

and 

hl s = (25) 

< V > t 

where < (< 77 > t ) < 00 is the time average of 77 = ^ over typical trajectory and 
i = 1, ...,n, n is the phase-space dimension. Transformation like u = ip{ x )i dr = r](x)dt 
is composed of a time re-parametrization followed by a space diffeomorphism. It is well 
known that the Lyapunov Exponents and Kolmogorov-Sinai Entropy are invariant under 
space diffeomorphism [T8] . 
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3 Lyapunov Exponent and Effective Potential: 



Now to compute the Lyapunov Exponent in terms of ISCO equation we shall first derive 
the radial potential in terms of Lyapunov Exponent. Therefore the Lagrangian for a test 
particle in the equatorial plane for any static spherically symmetric space-time can be 
written as 



<"5 



gtt i 2 + g rr r 2 + gcfxf, <fi 2 



(26) 



Now we defining the canonical momenta as 



* = f ■ (27) 

Using it, the generalized momenta can be derived as 

Pt = 9tti = -E = Const . (28) 
P<t> = g<t>4,<fi = L = Const . (29) 
Pr = g rr r . (30) 

Here (t, r, <p) denotes differentiation with respect to proper time(r). Now from the 
Euler-Lagrange equations of motion 

dp g _ SC 

Using it we get the non-linear differential equation in 2-dimensional phase space with 
phase space variables X^t) = (p r , r). 

dp r _ SC 

and 

£ = - ■ (33) 
dr g„ 

Now linearizing the equations of motion about circular orbits of constant r, we get the 
infinitesimal evolution matrix as 

Ma = ( I * ^ ) \r=r ■ (34) 

V grr J 

Now for circular orbits of constant r = r the characteristic values of the matrix gives 
the information about stability of the orbits. The eigen values of this matrix are principal 
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Lyapunov Exponent. Therefore the eigen values of the evolution matrix along circular 
orbits can be written as 



x 2 = J_! i 

g„ dr \5r J r r ° ' 

Again from Lagrange's equation of motion 

d_ /5£\ _ 6C _ 
d\ \Sr J Sr 



(35) 



(36) 



Therefore the Lyapunov Exponent (which is the inverse of the instability time scale asso- 
ciated with the geodesic motions) in terms of the square of the radial velocity (f 2 ) can be 
written as 



5C 



1 d 



(rg r 



Sr 2g rr dr 

Finally the Principal Lyapunov Exponent can be rewritten as 



A 2 lid 

2 g rr dr 



Again for circular geodesies [T7] 



r 2 = (r 2 )' = . 



(37) 



(38) 



(39) 



where prime denotes for a derivative with respect to r. Therefore the equation (1143ft must 
be reduced to 



or 



A 



A = ± 



'j.2\n 



(40) 



(41) 



We shall show in the next section, in any spherically symmetric spacetime the radial 
equation can be written as in terms of the standard effective potential which is given by 



r 2 = E 2 -V eff . 



(42) 



Therefore the Principal Lyapunov Exponent can be expressed in terms of standard effec- 
tive potential is 



A = ± 



V" 



eff 



(43) 



The Lyapunov Exponent must come in ± pairs to conserve the volume of phase space. 
Here VLt > determines the stability and V'Lf < indicates inst ability [9]. Alternatively 
the circular orbit is unstable and chaotic when the A is real, the circular orbit is stable 
when the A is imaginary and the circular orbit is marginally stable when A = 0. 

The above expression for A is valid for any stationary, spherically symmetric 
blackhole spacetimes i.e(Schwarzschild, Reissner Nordstr0m, Schwarzschild-de Sitter, 
Schwarzschild-Anti-de Sitter, Reissner Nordstr0m-de Sitter, Reissner Nordstr0m-Anti de 
Sitter etc.). 

4 Kolmogorov- Sinai Entropy and Lyapunov Expo- 
nent: 

An another important quantity which is related to the Lyapunov Exponents is so called 
Kolmogorov-Sinaij3] Entropy (hk s ), gives a measure of the amount of information lost or 
gained by a chaotic orbit as it evolves. Alternatively it determines how a system is chaotic 
or disorder when hk s > and non-chaotic for hk s = [18J. 

To compute it here, we use the formula which is given by PesinJI] that it is equal to 
the sum of the positive Lyapunov Exponents i.e 



In 2-dimensional phase-space, there are two Lyapunov Exponent, since h^s is equal to 
the sum of positive Lyapunov Exponent, therefore here the Kolmogorov-Sinai entropy in 
terms of effective potential is given by 



This entropy have played a crucial role in mathematical theory of chaos to check 
whether a trajectory in dynamical system is in disorder or not when it evolves with time. 
It is some sense different from the physical entropy, for example the entropy of the 2nd 
law of thermodynamics or blackhole entropy. Formally it is defined somewhat like entropy 
in statistical mechanics i.e it involves a partition of phase space. 

4.1 Critical Exponent and Effective Potential: 

Following Pretorius and Khurana|12j. we can define a Critical exponent which is the ratio 
of Lyapunov time scale T\ and Orbital time scale Tq can be written as 




(44) 



A;>0 




(45) 
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il T\ Lyapunov Time scale 
2tt\ Tq Orbital Timescale 

where we have introduced T\ = j and T w = ^f, which is important for black- hole merger 
in the ring down radiation. In terms of the square of the radial velocity (r 2 ), Critical 
Exponent can be written as 



T n 2vry (r 2 )" 

Now this exponent can be written as in terms of the inverse of the effective potential 
is given by 



^ = % = ^~w t ■ (48) 

Alternatively the reciprocal of critical exponent is proportional to the effective potential 
which is given by 



1 ^ _ _ V'lff 



^ = 2tt\ (49) 

5 Equatorial Circular Geodesies in Reissner Nord- 
str0m Spacetime: 

First, we shall consider a static, spherically symmetric, asymptotically flat solution of 
the coupled Einstein-Maxwell equations which represented by the Reissner Nordstr0m 
spacetime, the metric for such space-time becomes 



ds 2 = -^rdt 2 + ^dr 2 + r 2 (d9 2 + sin 2 6d<j) 2 ) . (50) 
r 2 A 

where A = r 2 — 2Mr + Q 2 . For M > Q, the horizons are situated at r± = M± \/ M 2 — Q 2 , 
for M = Q, the horizon is at r = M and for M < Q, there is no horizon, the singularity 
at r = is naked singularity 

To compute the geodesies in the equatorial plane for the Reissner Nordstr0m space- 
time, we follow ([17], [7J). To determine the geodesic motions of a test particle in this 
plane we set 9 = and 9 = constant = ~. 
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Therefore the necessary Lagrangian for this motion is given by 

1 



C 

2 



^ J-2 i r -2 , 2 12 

-— t + — r z + r z Z 
A 



(51) 



Now we defining the canonical momenta as 



5C . . 

P g = jr. (52) 



Using it, the generalized momenta can be derived as 



p t = ~t=-E = Const . (53) 



= r (f) = L = Const . (54) 

Pr = -r-r . (55) 

Here over dot denotes differentiation with respect to proper time(r). Since the Lagrangian 
does not depends on 't' and '0', so p t and are conserved quantities. Solving ( 1531) and 
(|54|) for t and we find 

i = E^. (56) 

= L/r 2 . (57) 

where E and L are the energy and angular momentum per unit mass of the test particle. 
Therefore the necessary Hamiltonian is given by 

H = p t i + p^cj) + p r r - C . (58) 

In terms of the metric the Hamiltonian is 

U = -^i 2 + T —r 2 + r 2 2 - C . (59) 
r A 

Since the Hamiltonian is independent of 't', therefore we can write it as 

2H = -^i 2 + r 4f 2 + r 2 <p 2 . (60) 
r z A 

= -Et + Lcj) + — r 2 = e = const . (61) 

Here e = — 1 for time-like geodesies, e = for light-like geodesies and e = +1 for space-like 
geodesies. Substituting the equations (1561) and (1571) in ( IBTl) . we obtain the radial equation 
for any spherically space-time is given by 

r* = i?-V e// = ^-(^- e U. (62) 
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where the standard effective potential for Reissner Nordstr0m spacetime is 

Veff = ~ e)^ 2 ■ (63) 



To investigate the circular geodesic motions of the test particle in the Einstein -Maxwell 
gravitational field , we must have for circular geodesies of constant r = tq and from the 
equation ( 1621) we get 



Veff = E 2 . (64) 

and 



dr 



. (65) 



Therefore we obtain the energy and angular momentum per unit mass of the test particle 
are given by 

= _ (rg-2Mr„ + Qy 

r 2 (rj - 3Mr„ + 2Q*) 1 ' 



and 



r 2 (Mr - Q 2 ) 



L r 2 -3Mr + 2Qi e ' (67) 

6 Non-Extremal Reissner Nordstr0m Spacetime 

6.1 Lyapunov Exponent and Equation of ISCO: 

For non-extremal case the radial potential of the test particle for time-like circular 
geodesies by substituting e = — 1 is given by 

L 2 \ ( 2M Q 2 



Veff = 1 + -7T 1 - — + 5r ■ (68 



and correspondingly the value of energy E and angular momentum associated with 
the circular orbit at r = tq are given by 

9 (r 2 - 2Mr + Q 2 f , s 

E = — f69l 

r 2 -3Mr + 2Q 2 ' 1 J 

2 r 2 (Mr - g 2 ) 

r 2 -3Mr + 2g 2 ' 1 } 
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Circular motion of the test particle to be exists when both energy and angular momentum 
are real and finite, therefore we must have r^ — 3Mr + 2Q 2 > and r > 

Also the angular frequency measured by an asymptotic observers for time-like circular 
geodesies at r = r is given by f2 

Q i v(Mr - m 

^0- j- -2 • \7l) 

Therefore for non-extremal Reissner Nordstr0m black-hole, the Lyapunov Exponents 
in terms of the ISCO equation becomes 



I -(Mrl - 6M 2 r 2 + QMQVq - 4Q^) 
Aw V l/ rg(r 2 - 3Mr + 2Q*) ' (?2) 

and the Kolmogorov-Sinai Entropy for non-extremal Reissner Nordstr0m black-hole be- 
comes 



-(Mrl ~ QM 2 r 2 + 9MQ 2 r - 4Q 4 ) 



r A (rl - 3Mr + 2Q 2 



So the time like circular geodesies of non-extremal Reissner Nordstr0m black-hole are 
stable when 

Mrl ~ §M 2 rl + 9MQ 2 r - 4Q 4 > . (74) 
such that \rn is imaginary, the circular geodesies are unstable when 



Mrl ~ §M 2 rl + 9MQ 2 r - 4Q 4 < . (75) 
i.e Xjin is real and the time-like circular geodesies is marginally stable when 

Mrl ~ §M 2 rl + 9MQ 2 r - 4Q 4 = . (76) 

such that Xrn is zero. 

The circular orbits are chaotic around the non-extremal Reissner Nordstr0m black-hole 
for Mr'l - 6M 2 r 2 + 9MQ 2 r - AQ 4 < i.e. h ks is real and for ISCO, h ks is zero. 

The solution of the preceding equation gives the radius of ISCO at r = rjsco for 
non-extremal Reissner-Nordstr0m black- hole [11] which is given by 

^ = 2 + Z 2 + ^. (77) 
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where 



Z 2 




1/3 



z 



1 



/ Q~ 2 Q* 



In the limit Q = 0, Z 3 = 4 and Z 2 = 2, we get the radius of ISCO, rjsco = 6M for 
Schwarzschild black-hole. 

Now the reciprocal of critical exponent in terms of ISCO equation for Reissner- 
Nordstr0m blackhole is given by 



For any unstable circular orbit, Tq > T\ i.e. Lyapunov time scale is shorter than the 
gravitational time scale making the instability observationally relevant [B]. 

6.1.1 Schwarzschild Black-hole: 

For Schwarzschild black hole a = Q = 0, the Lyapunov Exponents in terms of ISCO 
equation are 




(78) 




(79) 



and the Kolmogorov- Sinai Entropy for Schwarzschild blackhole becomes 




(80) 



So the time like circular geodesies of Schwarzschild black-hole are stable when 



r - 6M > . 



(81) 



such that Xsch or hk s is imaginary, the circular geodesies are unstable when 



3M <r <QM . 



(82) 
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i.e Xsch or or hk. s is real and the time like circular geodesies is marginally stable when 



r - 6M = . (83) 

such that Xsch or hk s is zero. This radius at r = r IS co — 6M gives the ISCO for 
Schwarzschild blackhole. The critical exponent for Schwarzschild blackhole is given by 



1 T n / -(r„ - 6M) 

7 = Tr 2 l (r„-3M) ' (84) 

For any unstable circular orbit say for ro = 4M, 7 becomes 2 J^_ , therefore Ta < Tq, 
i.e. Lyapunov time scale is shorter than the gravitational time scale, for this case the 
chaos will be damped and it is unobservable[6]. 

6.2 Lyapunov Exponent and Null Circular Geodesies: 

The radial potential that governs the null geodesies, can be expressed as by substituting 
e = in (|62D 



^-^-^-M). (85) 



where the effective potential for null geodesies is given by 

= ~ ■ (86) 
For circular null geodesies at r = r c , we know the condition from (|39|) 

U eff = E 2 . (87) 

and 

dUeff 



dr 



. 



Therefore we obtain the ratio of energy and angular momentum of the test particle eval- 
uated at r = r c for circular null geodesies are given by 



E = M-2«r. + q» (a) 



c 



r\ - 3Mr c + 2Q 2 = . (90) 
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After introducing the impact parameter D c = j?, the above equations reduced to 



1 



E 

T 



Mr c - Q 2 



(91) 



Solving equation (1901) we obtain the radius of null circular orbits are 

3M 



(92) 



It can be easily seen that from equation fl90|) circular null geodesies is the limiting case of 
the timelike circular geodesies. Therefore for this case the 2nd derivative of the effective 
potential is 



L 2 ( Q 2 



(93) 



Again the angular frequency measured by an asymptotic observers which is given by 

(94) 



n c = t = *n c = ±- 

t r* D c 



'Mr c - Q 2 



Using equation (19T|) we show that the angular frequency Q c of the circular null geodesies 
is inverse of the impact parameter D c . 

Using (|4T|) the Lyapunov Exponent for null circular geodesies is given by 



A 



RNN 




L 2 l 1 _2Q* 



(95) 



and the Kolmogorov- Sinai Entropy for Reissner Nordstr0m blackhole becomes 



h 



ks 




L- / 1 2Q 2 

'■; v r c 



(96) 



So the circular geodesies r c = r c+ and r c = r c _ are unstable since Xrnn is real. 
6.2.1 Schwarzschild Black-hole: 

For Schwarzschild black hole a = Q = 0, and the Lyapunov Exponents are 



A 



SchN 




(97) 



It can be easily check that for r c = 3M, Xs c hN is real, so the Schwarzschild photon sphere 
is unstable. 
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7 Extremal Reissner Nordstr0m Spacetime 

7.1 Lyapunov Exponent and Equation of ISCO: 

To derive the Lyapunov Exponent in terms of the equation of ISCO, we shall first write 
the effective potential for timelike circular geodesies which is given by 



r / V r J 



V eff = ( 1 + -T ) ( 1 " — ) • ( 98 ) 



Therefore we obtain readily the value of energy E associated with the circular orbit is 
given by 

2 (r - Mf 

r 2 (r - 2M) " ^ 

and the value of angular momentum associated with the circular orbit is given by 

L = . 100 

(r -2M) K J 

For circular motions, it must be required that energy E and angular momentum L, which 
are both real and finite, so we must have r — 2M > 0. Again the angular frequency 
measured by an asymptotic observers for time-like circular geodesies is given by Q 

t r 

Therefore for extremal Reissner Nordstr0m black-hole, the Lyapunov Exponents in terms 
of the ISCO equation becomes 

W = fp^ ^W^M)- (102) 

and the Kolmogorov-Sinai Entropy for extremal Reissner Nordstr0m black-hole be- 
comes 

ht , = 4^^y^lM } . (10 3) 
r oV( r o - 2M ) 

So the time like circular geodesies of extremal Reissner Nordstr0m black-hole are stable 
when 

r - 4M > . (104) 
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such that \ ex RN is imaginary, the circular geodesies are unstable when 



2M < r < 4M . (105) 
i.e \ ex RN is real and the time-like circular geodesies is marginally stable when 

r - AM = . (106) 

such that X ex RN is zero. 

Again the calculation for Kolmogorov- Sinai Entropy hk s for extremal Reissner Nord- 
str0m black-hole is valid for tq ^ r c . For r$ = r c = M, the hk s becomes zero. 

7.2 Lyapunov Exponent and Null Circular Geodesies: 

The effective potential that governs the null circular geodesies, can be expressed as by 
substituting e = in (1521) 

«.// = ^(i-f> 2 - (107) 

Therefore the ratio of energy and angular momentum of the test particle evaluated at 
r = r c for circular null geodesies are given by 



(108) 



E = ± (r c - M) 
L r 

r 2 c - 3Mr c + 2M 2 = . (109) 
After introducing the impact parameter D c = j?, the above equations reduced to 



Solving equation fl 109j) we obtain the radius of null circular orbits are 

r c = M, 2M . (Ill) 

It can be easily seen that from equation (11111) circular null geodesies is the limiting case of 
the timelike circular geodesies. Therefore for this case the 2nd derivative of the effective 
potential is 
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Again the angular frequency measured by an asymptotic observers which is given by 



t r% D c y 

Using equation (IllOp we show that the angular frequency fl c of the circular null geodesies 
is inverse of the impact parameter D c . 

Using (j4ip the Lyapunov Exponent for null circular geodesies is given by 




KxRNN = W-T 1 5- • (H4) 



So the circular geodesies r c = 2M are unstable since \ ex RNN is real. 

Note that for extremal blackhole, this result is valid for only single null geodesies i.e. 
r 7^ r c . For r = r c = M, the Lyapunov Exponent becomes zero i.e. \ ex RN = KxRnn = 0. 

Now we shall make a link between geodesic deviation equation and Lyapunov Expo- 
nent. 



8 Geodesic Deviation Equation and Lyapunov Expo- 
nent 

Geodesic deviation equation gives a proper measure of spacetime curvature. It connects 
the acceleration of the separation vector x 11 between two nearly geodesies. The separation 
four vector x{ T ) connects a point x a {r) on one geodesic(the fiducial geodesic) to a point 
x a (r) + x a { T ) on a nearby geodesic at the same proper time. In this section we shall find 
the link between Geodesic Deviation Equation and Lyapunov Exponent. 



8.1 Timelike Case: 

We know the timelike geodesic equation is given by 

(Px^ dx a dx@ 



1*7 = (115) 

where r is proper time along the geodesies and is the velocity four vector. Now 
to evaluate the geodesic deviation equation, x a {r) + x°( r ) must be obey the geodesic 
equation f)115p . Therefore we obtain the geodesic deviation equation which is given by 

d V , ore dx ° dxP I "tm dx a dx\ s 
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where \ 6 is the deviation four vector. We know the spherically symmetric RN spacetime 
can be written as 

ds 2 = -g{r)dt 2 + j^ dr2 + r 2 ( d ^ + sin 2 6d<p 2 ) . (117) 

where g(r) = 1 - 2M/r + Q 2 /r 2 . 

Using f)115p and (JTT7J) we get the following equations for timelike geodesies: 

dH = _g\r)_cttdr_ 
dr 2 g(r) dr dr 

1^ = -\ 9{r)9 ' {r){ t f + ^ (119) 
d 2 6 2drd9 „ „,d<b. 9 

d 2 2 dr d(f) ^ ^gd8d(f) (121) 

dr 2 r dr dr dr dr 

Hence for circular geodesies of constant r = ro on the equatorial plane 6 = n/2, these 
equations reduces to 



dH 

dr 2 



(122) 



>o)(|) 2 -<) 2 - (123) 

d 2 ^ 



dr 2 
d 2 </> 



(124) 

dT2 (125) 
and, using (I116p and (11171) we obtain following equations for geodesic deviation equation 

d 2 v* q'(ro) , dt w dv r N 

+ = < 126 > 

[5{(9'(ro)) 3 + 9(ro)9"(r„)}(^) 2 -{ 9 (r„)+r 9 '(ro)}(^) !! lx' = (127) 

£ + <g)V = (128) 
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dV 2 dd> ,d\ r , . 

-^ + ^ )( <JF> = < 129 > 

Again for timelike circular geodesies u.u = — 1 Therefore using d 1 1 T[) we get 

= (130) 

Solving equations (I123p and (I130p we have 

(131) 



,2 #'( r o) 



dr r[2g(r ) - r g'(r )} 

= ^ 

W 2g{r )-rg'{r ) 1 ; 

Using (11311) and (11261 - (11291) ) we transform the equation into (r — <p) plane we get the 
following equations 

def) 2 g{ro) dtp def) 
[\{{g'{r Q )f + g( ro )g"(r )}A 2 - {g(r ) + r g\r )}] X r = (134) 



d 2 Y e 

% + x' = (135) 



dV 2 dv r , . 

= < 136 > 

From the equation ( 11351) . we see that it represents a S.H.M, this implies that the 
motion in the equatorial plane 9 = it/ 2 is stable. Since the motion is harmonic, so we 
assume that the trial solutions of the remaining equations are 

X* = cie^ (137) 
X r = c 2 e^ (138) 
X* = c 3 e^ (139) 
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we get the reduced radial deviation equation is given by 

TV + 1 7W) 1 x = (140) 

Again for any static, spherically symmetric metric like (1117p the radial equation that 
governs the radial motions of the test particle is given by 

r 2 = E 2 - V eff = E 2 -(^-^ g(r) . (141) 

where the standard effective potential for any spherically symmetric spacetime 
(Schwarzschild, Reissner Nordstr0m) is 

V.// = " e ) • (142) 
and e has been previously defined. Once again we know for circular geodesies say r = r 

and 



Veff = E 2 . (143) 



dr 



. (144) 



Therefore we obtain the energy and angular momentum per unit mass of the test particle 
for timelike orbit are given by 

E 2 = - 2 f (ro)2 - - . (145) 
2g(r ) - rg'(r ) 

and 

T 2 r o9' '(r ) 



2g(r ) - r g'(r ) 



(146) 



Therefore the 2nd derivative of effective potential (j 142ft for timelike circular geodesies 
= — 1) of constant r = tq is given by 



V" ff 



9(r )g'(r ) + MnMrol _ 2(g'(r )f 



2g(r ) - r g'(r ) 



(147) 
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and the Lyapunov Exponent in terms of effective potential for any spherically symmetry 
spacetime is given by 



= V"iK (M8) 
Using (11401) and (I147j) we get the following equation in terms of the effective potential 



d 2 X r 



+ 



ro{2g(r ) - r g'(r )) 



VI ff X r = (149) 



2g'(r ) 

Using (I140p . (I147P and (JUS} we get the radial geodesic deviation equation in terms of 
the Lyapunov Exponent is given by 

d2 X r r (2g(r ) - r g'{r )) ^ 2 _ 

W %j sphX ~ ( } 

where the Lyapunov Exponent for any spherically symmetry blackhole is given by 



■ [g(ro)g'(r ) + 3 -^^-2(g'(r W] 
Kph ~ V 2g(r )-r g'(r ) " (151) 

Since 2g(r ) — r g'(r Q ) > and g'(r ) > for the allowed geodesic motion in the equatorial 
plane, so the solution of the equation is unstable for \ 2 sph > 0, stable \ 2 sph < and marginal 
stable for X 2 sph = 0. 

For Reissner Nordstr0m blackhole g(r ) = 1 — 2M/r + Q 2 /r 2 , 2g(r ) — r g'(r ) = 
r 2 — 3Mr + 2Q 2 and g'(r ) = 2(Mr ~Q 2 )/r^, therefore the Lyapunov Exponent becomes 



/ -(Mr 3 - QM 2 r 2 + 9MQ 2 r - 4Q^) 

which is similar to the expression ( 1721) and correspondingly the geodesic deviation equation 
is 

d 2 X T , r (Mr 3 - 6M 2 r 2 + 9MQ 2 r - 4Q 4 ) r 

2(Mr - Q 2 ) X ~ U 

For Schwarzschild blackhole g(r ) = 1 — 2M/r , 2g(r ) — rQg'(r ) = Tq — 3Mr and 
g'{ r o) = 2M/rQ, therefore the Lyapunov Exponent becomes 



which is similar to the expression ( 183]) and correspondingly the geodesic deviation equation 



is 



^ + ^_m xr = (155) 
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8.2 Null Case: 

For null circular orbit in the equatorial plane, the geodesic deviation is given by 



<i 2 Y M „ dx a dyP „ dx a dx 13 x 

Analogously for null circular orbit in the equatorial plane, the geodesic deviation equation 
is given by 

+ \[rl9"{r c ) - 2g(r c )} X r = (157) 
Again from (11421) the effective potential for null geodesies is given by 

U eff = £g(r) . (158) 
For circular null geodesies at r = r c 

Ueft = E 2 . (159) 

and 



dr 



. (160) 



Therefore we obtain the ratio of energy and angular momentum of the test particle eval- 
uated at r = r c for circular null geodesies are given by 



E g(r c ) 

T = ± \—- (161) 



2g{r c ) - r c g'{r c ) = . (162) 

It can be easily seen that from equation (I162p circular null geodesies is the approxi- 
mation of the timelike circular geodesies. 

Therefore for this case the 2nd derivative of the effective potential is 

K ff = ^(r 2 c g"(r c -2g(r c )) . (163) 

' c 

Thus the Lyapunov Exponent in terms of effective potential for null circular geodesies 
can be written as 



vv„// = \/ • (164) 
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Using (I157p and (11471) we get the following equation in terms of the radial potential 
becomes 

d 2 y r ( r 4 \ 

+ U7i^// \r=r<X r = (165) 



d<\> 2 ^ \2L 2 eff J {r=r 

Using (11571) . ( I147p and (11641) we get the radial geodesic deviation equation in terms of 
the Lyapunov Exponent for null case is given by 

rfV fr f 2 >?Nuu)x r = (166) 
where the Lyapunov Exponent for null circular geodesies is given by 



\Null 



■^[2g(r c )-r^g"(r c )] . (167) 



The solution of the equation is unstable for \ 2 Nuil > 0, stable \ 2 Nu u < and marginal 
stable for = for the allowed geodesic motions in the equatorial plane. 

For Reissner Nordstr0m blackhole g(r c ) = 1 — 2M/r c + Q 2 /r 2 and 2g(r c ) — r 2 g"(r c ) = 
2(1 — 2%) therefore the Lyapunov Exponent becomes 



*Nuii = \/^(l-2^) . (168) 

which is similar to the expression ( I95j) and correspondingly the geodesic deviation equation 
is 

^ U-2^j X r = (169) 

For Schwarzschild blackhole g(r c ) = 1 — 2M/r , 2g(r c ) — r 2 g"(r c ) = 2 therefore the 
Lyapunov Exponent becomes 




Xnuii = \ — ■ ( 17 °) 
y ' c 

which is similar to the expression ( 1971) and correspondingly the geodesic deviation equation 
is 



-,2 



-X r = (171) 
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9 Ratio of Angular velocity of time like circular orbit 
to Null Circular Orbit 



Since we have already proved that for time-like circular geodesies the angular velocity is 
given by from equation (|71|) 



_ y/(Mr - gj 
"0 — 5 



(172) 



Again we obtained for circular null geodesies fl c = 4-, so we can deduce similar expression 
for it is given by 

' c 

Resultantly we obtain the ratio of angular frequency for time-like circular geodesies to 
the angular frequency for null circular geodesies is 




(174) 



For tq = r c , Qq — ^ C) i- e > when the radius of time-like circular geodesies is equal to the 
radius of null circular geodesies, the angular frequency corresponds to that geodesic are 
equal, which demands that the intriguing physical phenomena could occur in the curved 
space-time, for example, possibility of exciting Quasi Normal Modes(QNM) by orbiting 
particles, possibly leading to instabilities of the curved space-time [7]. 

For r > r c , we shall prove that for Schwarzschild black-hole, Reissner Nordstr0m 
black-hole the null circular geodesies have the largest angular frequency as measured by 
asymptotic observers than the time-like circular geodesies. We therefore conclude that 
null circular geodesies provide the fastest way to circle black holes[8j. 

Now the ratio of time period for time-like circular geodesies to the time period for null 
circular geodesies is given by 



T / y/Mr c - Q 



T c V J Mr, - Q2 




(175) 



This ratio is valid for ro ^ r c . For ro = r c , To = T c , i.e. time period of both geodesies 
are equal, which possibly leading to the excitations of Quasi Normal Modes. For r > r c , 
T > T c , which implies that the orbital period for time-like circular geodesies is greater 
than the orbital period for null circular geodesies. For r = r IS co an d r c = r photon , 
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therefore the ratio of time period for ISCO (r = risco) to the time period for photon- 
sphere (r c = r photon ) for Reissner Nordstr0m black-hole is given by 



Tisco = (V^EV)(d) . (176) 
T phot0 n \y/Mr -Q 2 J VI J 

This implies that T ISC o > T photon , therefore we conclude that ISCO provide the slowest 
way to circle the Reissner Nordstr0m black-hole. 



9.1 Schawrzschild Black hole 

For Schawrzschild black hole, a = 0, Q = we get the ratio of angular frequency for 
time-like circular geodesies to the angular frequency for null circular geodesies is 



(177) 



which is proportional to the radial coordinates r. Since for Schwarzschild black hole, the 
ISCO occurs at rjsco > 6M and photon sphere occurs at r c = 3M, therefore the ratio of 
angular frequency for ISCO(r = r IS co — 6M) to the angular frequency for photon-sphere 
(r c = r phot0 n = 3M) is given by 



n isco = r c ^F c = J_ ^ 

^photon 

r 0v /r^ 2\/2 ' 

Which suggests that, Qisco < ^photon, i-e ISCO is characterized by the smallest angular 
frequency as measured by the asymptotic observers than the null circular geodesies. Cor- 
respondingly the ratio of time period for time-like circular geodesies to the time period 
for null circular geodesies is given by 

£ = (179) 



Since for Schwarzschild black hole, the ISCO occurs at rjsco ^ 6M and photon sphere 
occurs at r c = 3M, therefore the ratio of time period for ISCO (r = r IS co — 6M) to the 
time period for photon-sphere (r c = r photon = 3M) is given by 

Tjsco _ W^ = 2v ^. (lg0) 



T, 



photon > cy > c 

This implies that Tjsco > Tph oton , i.e the orbital period for time-like circular geodesies 
is greater than the orbital period for null circular geodesies. From this we conclude that 
innermost stable circular orbit provide the slowest way to circle the black hole. 
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10 Discussion 



In this article, we have used Lyapunov Exponent and Kolmogorov- Sinai entropy to give 
a full description of time-like circular geodesies and null circular geodesies in Reissner 
Nordstr0m black hole space-time. We then explicitly derived them in terms of the equation 
of the radius of ISCO. We then further showed that the Lyapunov Exponent can be 
used to determine the stability and instability of equatorial circular geodesies, both time- 
like and null case for Reissner Nordstr0m black hole space-time. Analogously, we also 
compute the Kolmogorov- Sinai entropy to measure the disorder-ness of a chaotic orbit as it 
evolves. Using them, we found that the ISCO occurs at risco — 4M for extremal Reissner 
Nordstr0m black-hole and rjsco — 6M for Schwarzschild black-hole. We also made a link 
between Lyapunov Exponent and Geodesic Deviation equation. The other point we have 
studied that for circular geodesies around the central black-hole, ISCO is characterized by 
the smallest angular frequency as measured by the asymptotic observers-no other circular 
geodesies can have a smallest angular frequency. Thus such types of space-times always 
have flparticie < ^photon for all time-like circular geodesies. Alternatively it was shown that 
the orbital period for time-like circular geodesies particularly, ISCO is characterized by 
the greatest orbital period than the null circular geodesies. Therefore, we conclude that 
Innermost Stable Circular Orbit (ISCO) provide the slowest way to circle the black hole. 
The aim in future to find the ISCO for extremal Kerr and Kerr-Newman by using the 
Lyapunov Exponent and Kolmogorov- Sinai entropy. 
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